Multibrane Inflation and Dynamical Flattening of the Inflaton Potential 
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We investigate the problem of fine tuning of the potential in the KLMT warped flux compactifi- 
cation scenario for brane-antibrane inflation in Type IIB string theory. We argue for the importance 
of an additional parameter ipo (approximated as zero by KLMT), namely the position of the an- 
tibrane, relative to the equilibrium position of the brane in the absence of the antibrane. We show 
that for a range of values of a particular combination of the Kahler modulus, warp factor, and ipo, 
the inflaton potential can be sufficiently flat. We point out a novel mechanism for dynamically 
achieving flatness within this part of parameter space: the presence of multiple mobile branes can 
lead to a potential which initially has a metastable local minimum, but gradually becomes flat as 
some of the branes tunnel out. Eventually the local minimum disappears and the remaining branes 
slowly roll together, with assisted inflation further enhancing the effective flatness of the potential. 
With the addition of Kahler and superpotential corrections, this mechanism can completely remove 
the fine tuning problem of brane inflation, within large regions of parameter space. The model can 
be falsified if future cosmic microwave background observations confirm the hint of a large running 
spectral index. 

PACS numbers: 11.25. Wx, 98.80.Cq 



I. INTRODUCTION 

The past few years have witnessed significant progress 
in the building of realistic inflation models within string 
theory 0. The most popular approach has used the in- 
teraction potential between branes and antibranes as the 
source of energy driving inflation; hence the interbrane 
separation plays the role of the inflaton. Early attempts 
suffered from the assumption that moduli could be fixed 
by fiat, without specifying the stabilization mechanism. 
However, explicit examples have shown that properly fix- 
ing the moduli can interfere with the flatness which is re- 
quired for the inflaton potential 2] (KLMT). Flux com- 
pactifications in type IIB string theory have provided a 
robust way of stabilizing the dilaton and complex struc- 
ture moduli of the Calabi-Yau manifold |3, with the 
added benefit of providing warped throats {Jf which can 
be useful for tuning the scale of inflation or of the stan- 
dard model. Nonperturbative effects have been plausibly 
invoked to stabilize the remaining Kahler moduli 0, |(| . 

This last step, unfortunately, generically spoils what- 
ever flatness of the inflaton potential which was accom- 
plished by warping, and thus requires additional fine tun- 
ing of parameters. One can introduce extra parameters 
associated with corrections to the superpotential, which 
must be tuned to a part in 100 in order to obtain the 
minimum amount of 60 e- foldings of inflation 2j. Alter- 
natively, it is possible to tune parameters already present 
within the model but at the level of a part in 1000 0. 

The fine-tuning problem of the KLMT model has been 
considered by several authors H, . Shift symmetry of 
the Lagrangian has been proposed as a possible solution 
|10|. It has also been pointed out that if Lyman a con- 
straints are viewed with skepticism, the degree of tuning 
required is considerably relaxed [ll|. Of course, this ob- 



servation also applies to other models of inflation. 

In this paper we reexamine the fine-tuning issue, start- 
ing within the context of the original KLMT potential, 
without introducing additional corrections. A closer con- 
sideration of the potential between the brane and an- 
tibrane reveals the importance of a parameter (which we 
call ipo ) which was set to a certain value in KLMT. It is 
the distance which the brane will travel before annihilat- 
ing with the antibrane in the warped throat, assuming 
that the brane starts somewhere closer to the the inte- 
rior of the Calabi-Yau at an initial position ip = Q. The 
other important parameters governing the shape of the 
inflaton potential are the warp factor of the throat, which 
we denote as e 1 / 4 , and the Kahler modulus at its stabi- 
lized value, a. The warp factor is tunable through the 
ratio of NS-NS and RR fluxes K,M: e 1 / 4 = e -^KI?>9«M 
0, 3 . The VEV of the Kahler modulus a is determined 
by the nonperturbative superpotential W — Ae~ aa in- 
voked by ref. |5j, where a = 2ir/N in the case of gluino 
condensation in an SU(N) gauge theory sector which can 
naturally arise in these compactifications. 

If the parameter ipo is set to zero as in KLMT, then 
indeed it is necessary to introduce additional sources of 
dependence on the inflaton into the potential, in order 
to cancel the unwanted contribution to the inflaton mass 
coming from the Kahler modulus. However, in this work 
we point out that for a range of values 0.259 < b < 0.286 
of the combination b 2 = 2cre?pQ 6 , it is possible to achieve 
60 or more e-foldings of inflation without any additional 
corrections to the potential. Fortunately this restriction 
on b applies only to the simplest KLMT model. We will 
show that the restriction disappears when one considers 
more general examples in which corrections to the Kahler 
and superpotentials are taken into account. 

Another crucial parameter which we consider is the 
number of mobile mobiles branes, N, generalizing the 
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simplest case of only a single brane playing the role of 
the inflaton. We refer to this possibility as multibrane in- 
flation. Within the aforementioned range of b values, we 
observe a very interesting phenomenon, whereby start- 
ing with a sufficiently large number N dynamically leads 
to a sufficiently flat inflaton potential. At large N, the 
branes start out being confined in a metastable minimum 
of the potential. But as successive branes tunnel out of 
this minimum, it becomes shallower, until finally all the 
remaining branes can roll together into the throat. This 
is illustrated in figure H The potential is often very flat 
when the branes start rolling, leading to a long period of 
inflation. This is a novel, dynamical mechanism for flat- 
tening the inflaton potential, which appears to be quite 
specific to brane-antibrane inflation. 
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FIG. 1: Sequence of inflaton potentials due to successive tun- 
neling of branes toward the right. 



II. WHEN ARE MULTIPLE BRANES 
ALLOWED? 

Before looking at the dynamics of inflation with stacks 
of branes and antibranes, we consider under what condi- 
tions multiple branes are compatible with the stabiliza- 
tion of the extra dimensions, i.e., the Kahler modulus, 
(j, in supergravity language. For each extra brane we 
have to introduce a corresponding antibrane such that 
the tadpole cancellation condition is statisfied. Adding 
antibranes lifts the local minimum of the potential by an 
additional amount, given by the sum of their tensions, 
reduced by the warp factor at their location in the bot- 
tom of the Klebanov-Strassler throat. There is a critical 
number of antibranes that can be added such that the 
potential still has a local minimum along the a direc- 
tion and therefore the Kahler modulus remains stable. 
The stronger the warping, the larger the number of an- 
tibranes that can be added. In our analysis we will add 



antibranes such that the Kahler modulus remains sta- 
bilized and heavy. Thus it does not act like a second 
inflaton field, since it changes very little during the slow- 
rolling of the brane-antibrane separation if). 

Not all models are suitable for constructing such a 
multi-brane setup. We will give here two examples, one 
which works and one which does not. The original KKLT 
setup does not accommodate a large number of mo- 
bile branes; just one additional antibrane is sufficient to 
lift the vacuum from anti-deSitter (AdS) to Minkowski 
(or deSitter, dS). The barrier that prevents tha Kahler 
modulus from rolling has a height roughly equal to the 
depth of the AdS minimum. Addind a second brane lifts 
the potential by an amount that is roughly equal to the 
height of the barrier, which destroys the local minimum. 
In figure|21we take the parameters of the potential and the 
antibrane tension as given by .S]. One sees that adding 
a second brane results in a very shallow minimum and 
three branes destabilizes the Kahler modulus. The po- 
tential of the Kahler modulus is given by 

n Af~ aa /i \ (VD 

V = ^a^\i e - aa + Wo+Ae- a °)+^ (1) 
with W = -f(T 4 , A = 1, a = 0.1 D = 3 X 1(T 9 . 




FIG. 2: Potential energy of the Kahler modulus in the the 
original KKLT scenario for N = 0,1,2 and 3 antibranes, 
showing there can be at most 2 antibranes without desta- 
bilizing the Kahler modulus. 

The situation is dramatically im prov ed if one uses a 
racetrack superpotential as in ref. |l2j| . The minimum 
is already Minkowski and adding antibranes in a highly 
warped throat lifts the vacuum to dS without destabiliz- 
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ing the Kahler modulus. The potential is: 

aAe hG ) 



-2(a+b)a 

V = — (bBe ao 



6<r 2 



(2) 
ND 



[Be aa {3 + ba) + e ba (A(3 + aa) + 3e aa W )] 



where we choose the same parameters as in 12], A — 
1, B = -1.03, a = 2tt/100, b = 2tt/99, W = -2 x 10~ 4 
and D — 10~ 9 . With these values the critical number of 
branes turns out to be N — 28, illustrated in fig. 01 



N = 29 



N =28 



FIG. 3: Using a racetrack superpotential and a throat with 
stronger warping one can add a large number (28 in this exam- 
ple) of antibranes while keeping the Kahler modulus stable. 

In the remaining analysis, we will not be concerned 
with the particular model of Kahler modulus stabiliza- 
tion, beyond the assumption that it admits the required 
number of additional branes. It will also be found that 
a large number is not required, so the racetrack model 
mentioned here is sufficient to establish the existence of 
a viable model. 



III. INFLATON ACTION 

Our starting point for the dynamics of inflation is the 
Lagrangian for the mobile brane position ip in the KLMT 
framework. To simplify the analysis, we will assume that 
all moduli (notably the Kahler modulus a) are heavy 
except for the inflaton, and that the 4D effective cos- 
mological constant vanishes at the end of inflation. The 
Lagrangian can be parametrized as 



r 6(7 12 Te 



where r — (2a — if) 2 ) ^jjj, T is related to the tension r of 
the 3-branc or antibrane by T = 2r/g 4 s @, e is the warp 
factor, ip is the 3-brane position, and ipo is the location 
of the antibrane. To keep our work self-contained, this 
expression is derived in Appendix A, in just the same 
way as in KLMT except for the omission of two approx- 
imations: (1) we do not assume that ip 3> ipo, and (2) we 
do not Taylor-expand the term (1 + e/(ip — i/jq) 4 )^ 1 ■ 

Let us first comment on point (1) above. In KLMT 
and certain other investigations of brane-antibrane in- 
flation, it was assumed that inflation takes place suffi- 
ciently far from the antibrane to justify the approxima- 
tion \ip — ipo\ — V'o- On the other hand, the ipo parameter 
was kept in Q and 0]. Physically, it may be impor- 
tant to allow for nonzero ipo because of the possibility 
that the ■i/'-dependence coming from the (2a — ip 2 ) 2 fac- 
tor in the potential will tend to stabilize the brane at 
some point other than the bottom of the warped throat, 
in the absence of the brane-antibrane interaction. An 
example where this is particularly clear is that of two 
symmetrically placed warped throats, each with its own 
antibrane [14(- In that case the potential (translated to 
our notation) was taken to be 

V^\2- 



[i> - Vo) 4 (V> + 



(4) 



1 + 



(3) 



From symmetry considerations it is obvious in this model 
that there exists a metastable or unstable equilibrium 
point exactly midway between the two throats, at ip = 0, 
which cannot coincide with the bottom of either throat. 
If we simply remove one of the throats, the result corre- 
sponds to our starting point 

The presence or absence of the t/'o parameter explains 
differences in the amount of required fine-tuning of the 
potential found in different investigations. To elucidate 
this, recall that the ^-dependence in the volume modulus 
r = (2a - ip 2 ) is the origin of the r\ problem pointed out 
in KLMT — the mass of the inflaton is of order W ~ H 
in Planck units, ruining the slow roll condition. It was 
noticed by [jj that it is possible to tune parameters of 
a model with a potential similar to that in © so that 
the slow-roll conditions are satisfied, even without adding 
any extra ^-dependent corrections to the superpotential, 
by varying only the parameters which already appear. 
However, the degree of tuning needed was more severe 
than in KLMT where corrections to the superpotential 
were invoked. One way of understanding this is to realize 
that if tpo — and only corrections of order ip 2 are added, 
then the potential automatically has V'(Q) — 0, so that 
one of the slow roll parameters V /V is small without 
any tuning. Then only the r\ parameter (V" /V) needs 
to be tuned. In our potential however, both slow-roll pa- 
rameters must be tuned, for generic parameters in the 
potential. In it was found that 60 e-foldings of infla- 
tion could be obtained by tuning the warp factor to one 
part in 2000, whereas KLMT only needed one part in 
100 tuning. Although this sounds contrary to the spirit 
of the present paper, it was due to the fact that no sys- 
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tematic study of the effect of varying -0o was done in ; 
instead only one fixed value was considered. 

We note that the behavior of the potential as ip —> 4>o 
is classically correct, since by derivation it is of the form 
l/h(ip) where h(ip) satisfies V 2 h — C8(ip — ipo) in the 6D 
Calabi-Yau manifold. Although the full brane position 
must be specified by 3 complex ip coordinates, we are 
keeping the angular ones fixed and following only the 
radial coordinate for the brane-antibrane separation. 

We also comment on our choice of avoiding the second 
approximation (point (2) above), by not Taylor expand- 
ing (1 + e/(ip — i/jo) 4 )~ X - We choose not to do so be- 
cause there is no advantage in making this expansion — it 
only leads to an unphysical divergence in the potential 
energy as the brane-antibrane distance shrinks. The un- 
expanded form has good behavior as ip — > ipo, including 
the vanishing of the potential at the minimum. 



A. Multibrane Lagrangian 

We will be interested below in the generalization to 
more than one brane rolling. The Lagrangian for N 
branes and N antibranes is 



„ 6cr A .„ NTe ( A e 



-1 



(5) 



where now the volume modulus is 

N 

r = 2a-Yj\ (6) 



Later we will specialize to the situation where all the 
branes are coincident, tpi = ip. Then the potential takes 
the form NeT r - 2 (l + Ne(ip - ipo)' 4 )' 1 . 



IV. CONDITIONS FOR INFLATION 

The optimal potential for achieving slow roll inflation 
is one where both the first and second derivatives vanish 
simultaneously at some critical point ip c . Although this 
requires fine tuning, we will discuss a mechanism where 
the tuning can automatically occur through a dynamical 
process. To characterize the tuning, we introduce new 
parameters 



_ Ne b2 _ 2<re _ ip 
such that the potential is proportional to 



V(x) 



1 



(b 2 1 a — x 2 y 



1 



(x-iy 



(7) 



(8) 



The conditions for the potential to have a flat point, as 
defined above, are 



These equations can be solved analytically, as the so- 
lution of a cubic equation, but the expression is not very 
enlightening. Instead we give an approximate descrip- 
tion, by first considering b as a free parameter, and de- 
termining as a function of b the critical values a c and 
x c which satisfy ©■ For consistency, we should also re- 
quire that b 2 > a; otherwise the potential blows up at 
x = b/^/a < 1, before the brane reaches the antibrane. 
The range of values where a c > (obviously also required 
for a physically sensible solution) and a c < b 2 turns out 
to be rather narrow, 



0.259 < b < 0.286, 



(10) 



and in this region a c is to a good approximation a linear 
function of 6: 



= 36 - 0.77 



(11) 



One can also approximate the critical value x c by a 
quadratic form, x c = -0.146 + 1.215 + 6(6 - 0.259) (b - 
0.286). 

We will show later that the narrowness of the inter- 
val i|10|) is an accident, in the sense that physically well- 
motivated perturbations to the potential can remove this 
restriction. Therefore we will ultimately not regard it as 
a fine tuning. However the special value a c for a given 
value of b clearly does represent a fine tuning of the pa- 
rameters. It is this condition which will be addressed 
by a dynamical mechanism. Namely, by starting with a 
large enough number N of branes, one can initially have 
a > a c (but still a < b~). The potential always has a 
metastable minimum from which branes can tunnel in 
this case. After each tunneling event, the value of a di- 
minishes by a discrete amount. At some critical number 
of branes, N c , the curvature changes so that the branes 
roll instead of being trapped. It becomes a quantitative 
question, depending on the values of a and e, whether 
the potential is sufficiently flat to give enough inflation. 
We will show that in fact it is almost always flat enough 
to give at least 60 e-foldings of inflation. 

To find out how likely it is that adequate inflation will 
occur, we have explored the parameter space e, a, b, 
subject to the aforementioned restrictions and the as- 
sumption that N starts out being larger than the critical 
value, where all the branes start to roll. For each value 
of {e, a, 6}, the quantities ipo an d N c are determined. 
The fields start from rest at the value where the poten- 
tial has a local minimum with N = N c + 1 branes, on the 
assumption that the nontunneling branes do not move 
during the tunneling event. The equations of motion are 



~dt 



dV 

dipt' 



dNe 
dt 



H 



(6x- l)(x- l) 4 , 5x 2 (x-l) 4 



(9) 



for the canonical momentum 

12cr 

7Tj = — 5- 



1>i 



(12) 



(13) 



and number of e- foldings N e = In a(t), where the Hubble 
parameter is given by H 2 = (6aNip 2 /r 2 + V)/(3M 2 ). 
Since all the fields move together, we need consider only 
one equation for the inflatons, which we do numerically. 

The result is that, for b satisfying 1)10(1. there is a high 
probability of getting enough inflation, for random values 
of e, a selected from within a reasonable range, —5 < 
log 10 e < —2.4, 1 < log 10 2a < 3.6. Scanning uniformly 
over this region, we find that 85% of tries result in more 
than 55 e-foldings of inflation when b = 0.27, and an 
even greater percentage occurs for other values of b. The 
results are illustrated for b = 0.27 in Figure 01 showing 
the number of e-foldings, N e (omitting values where N e > 
300 for clarity) versus log 10 e and log 10 2a. 

In figure^it can be seen that N e is constant along lines 
of constant A = (2cr) 2 / 3 e _1 ' 3 . This can be understood 
by rewriting the potential in terms of the canonically 
normalized field cj> = tpiy/3N/a: 



N, 



v 



TN 



1 



1 - 



36(N/Af 

(<t>-4>o) 4 



(14) 



where <po — b~ x l 3 y/3N/A. Thus for fixed &, the shape of 
the potential, power spectrum, and duration of inflation 
depend only upon combination of parameters in A, and 
not on the orthogonal combination B = (2cr) 1 / 3 e 2 / 3 . We 
can therefore restrict our exploration of parameter space 
to a one-dimensional scan along the direction of A. More- 
over, we can show that the number of branes TV depends 
in a very simple way on the value of A: N — ab^ 4 / 3 A. 
This follows from eliminating ipo from the definitions of 
a and b. Therefore the critical number of branes is given 
by 



N r = a r b- 4 / 3 



(2a) 2 / 3 
,1/3 



(15) 



and for a fixed value of 6, we can treat the parameters A 
and N c interchangeably. Of course H15|) can't be satisfied 
exactly, because N c has to be an integer. The acciden- 
tal discrepancies between the optimal value of N c and 
the closest integer below this value will affect how long 
inflation actually lasts. 

In figure |3] we show the correlation between number 
of e-foldings and the critical number of branes, scanning 
uniformly in the interval 1 < log 10 A < 4 of parameter 
space. Although the two quantities are generally corre- 
lated, there are order-of-magnitude variations in N e over 
narrow ranges of N c . Figure |SJ demonstrates that both 
features can be matched fairly well by an empirical rela- 
tionship. Let AN C be the discrepancy between the ideal 
value of N c and the closest integer below this value. Then 
we find that 



iV e = 10 



AN C 



(16) 



The factor l/AN c , is understandable since if AN C = 0, 
the potential is tuned to be perfectly flat at some position 
tp = tp c . 
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FIG. 4: iV e versus log 10 e and log 10 2a for b = 0.27, in the 
range iVe < 300. 
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FIG. 5: log iVe versus the critical number of branes log N c 
b = 0.27. The empirical relationship 1161 is also plotted. 



The factor \/W~ c in H16|l can be partly understood in 
terms of assisted inflation ^{|. With N fields rolling si- 
multaneously, the kinetic term is multiplied by N relative 
to one field. We should rescale the fields by tp — ► rp/ y/N, 
and the slow-roll parameters by 1/N. One therefore ex- 
pects an enhancement in the length of inflation by some 
power of N c . It would of course be satisfying to under- 
stand why that power is 1/2 in the present model. Wc 
have not yet succeeded in analytically deriving this result. 
However it is not surprising that N e fails to be enhanced 
by a full factor of N c as would be the case for ordinary 
assisted inflation, because the potential contains Nip 2 in 
the Kahler modulus. Only if the potential was indepen- 
dent of N would we be sure that N e ~ N c . In that case, 
transforming to the canonically normalized field would 
cause V(ip) — > V(<f>/y/~N) and the slow roll parameters 
would go like 77 ~ V" /V ~ 1/N. The surprise is that 
there is any assisted inflation effect at all. The reason for 
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FIG. 6: log N e versus log iV c for b = 0.26, 0.27, and 0.28. 
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FIG. 7: Log of the inflationary scale as a function of the 
number of branes. 



it is that at the critical value of 7Y, the potential is flatter 
than it would generically be for an arbitrary function of 
Nip 2 . 

The above results were for b = 0.27. Allowing b to vary, 
ones finds that smaller b tends to give more inflation with 
a smaller number of branes. This dependence is shown 
in figure |H| The main observation, though, is that the 
great majority of all cases yield more than 60 e-foldings 
of inflation. 



V. SCALAR POWER SPECTRUM 

In addition to having a long enough period of inflation, 
we must obtain a spectrum of density perturbations in 
agreement with observational constraints from the Cos- 
mic Microwave Background. To compute the power spec- 
trum, we use the following expression |l6j ] which applies 
to a set of noncanonically normalized fields with La- 
jfripi - V: 



jrangian C = | Gy ( 



P(k) 



V dN e dN e 
75tt 2 M 2 dip 1 dt/ji 



Qij 



N C VH 2 r 
75ir 2 M 2 6^2 



(17) 



These quantities are evaluated at horizon crossing of the 
relevant wave number, when k = aH . To normalize the 
spectrum, we evaluate P(k) for modes which crossed the 
horizon 55 e-foldings before inflation ended, correspond- 
ing to a scale of inflation around I/ 1 / 4 ~ 4 x 10 14 GeV. 
The consistency of this assumption will be verified. The 
COBE normalization implies that P = 4x 10 -10 for these 
modes. This allows us to normalize the scale of the po- 
tential by fixing the value of the parameter T which is 
proportional to the unwarped brane tension. The results 
are expressed in terms of the scale of inflation relative 
to the Planck scale. In figure we plot log(V 1/4 /M p ) 
as a function of the number of branes. (In the previous 
section it was explained that, just like N e , spectral prop- 
erties depend only on the same combination of a and 



0.15 



0.1 



0.05 



-0.05 



b 


= 0.28 


o b 


= 0.27 


= b 


= 0.26 





J <pWJW ,1 WV>^o o «»o o °o O0o0 <>oC> o ° o o o o o o o a 



N, 



400 



600 



FIG. 8: Deviation of scalar spectral index from unity versus 



e which determine N c .) From the figure one sees that 
the inflationary scale ranges between 10~ 3 and 10 Af p , 
roughly consistent with our assumption of ~ 55 e-foldings 
of inflation. The string scale, which determines the brane 
tension, must be in about the same range. If the number 
of branes is small enough, the inflationary scale can be 
sufficiently high for the tensor contribution to the cosmic 
microwave background to be observable in future exper- 
iments. 

The scalar spectral index is given by n s = 1 + 
d\nP/d\nk. We have numerically evaluated it over the 
relevant range of parameters. Like the other quantities, it 
is most conveniently displayed as a function of the num- 
ber of branes. The deviation n s — 1 at 55 e-foldings be- 
fore the end of inflation is shown as a function of N c 
in figure |SJ There is a wide range of possible values, 
-0.07 < n s - 1 < 0.15. 

It is also interesting to consider the running of the 
spectral index, dn/d\nk = d 2 In P/ dink 2 , for which 
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FIG. 9: dn/dlnk versus n s 



FIG. 10: dn/d In fcwithn-1, with WMAP+SDSS constraints 
shown. 



WMAP has given a hint of a larger-than-expected neg- 
ative value, dn/d\nk = —0.02 |17|. Although the er- 
ror bars are large enough to be consistent with no run- 
ning, the suggestion of large running has attracted at- 
tention. Generically, within the slow-roll formalism, one 
expects dn/dhik ~ (n s — l) 2 , which is usually smaller 
than the measured central value. Our numerical survey 
shows that it is possible to have rare instances of run- 
ning as negative as —0.01, even though inflation lasts 
long enough. However, as shown in figure El it is still 
true that dn/d\nk ~ (n s — l) 2 . The rare cases of large 
running correspond to occasional extreme values of the 
spectral index. 

We have checked that the spectral index and running 
are in good agreement with the predictions of slow roll 
inflation, n—l = 2r] — 6s, dn/d\nk = — 2£+ 16erj — 24e 2 , 
where e = \(V /V) 2 , rj = V" /V and f = V'V"'/V 2 , 
in units where M p = 1. The derivatives must be taken 
with respect to the canonically normalized fields, <pi = 
(r/\/12a)ipi. In all cases, e ~ 10~ 6 — 10~ 9 , whereas 
i] ~ —0.03 and £ ~ 10~ 3 . The running is completely 
dominated by £. 

It is interesting to compare the predicted values of the 
running and the spectral index with experimental con- 
straints from combined WMAP and Sloan Digital Sky 
Survey Lyman a data |l8j The allowed region is super- 
imposed on the predictions in figure ITU1 Although most 
of the predicted points are allowed, we note that the few 
examples with large running have too large a deviation 
of n s from unity. Thus if large running should be ex- 
perimentally confirmed in the future, it will rule out the 
models under consideration. 



VI. GENERALIZED MODELS 

The previous analysis shows that within a narrow 
range of the parameter b = y/2ae/ipQ, a long period of 
inflation will always ensue, as long as the initial number 



of mobile branes exceeded the critical number N c . But 
since b looks fine-tuned, this is not a complete solution 
to the problem of naturalness for the inflaton potential. 
Could this be an accident of the particular form of the 
simplest KLMT potential? To answer this, we explore 
the effect of perturbing the model by corrections which 
are expected to be present, but which we have so far 
ignored for the sake of simplicity. 

At least two kinds of modifications to the basic KLMT 
model are generically expected: corrections to the super- 
potential and to the Kahler potential. A simple example 
of the former was considered in Q , where the superpo- 
tential for the Kahler modulus T = a + iO and inflaton 
was taken to be 

W — W + g(T)(l + Sip 2 ) (18) 

leading to a contribution to the potential of the form 

a + f3i! 2 + 7V> 4 



sv 



(19) 



with CT-dependent coefficients a, (3, 7. This should be 
added to the brane-antibrane potential which we have 
been using so far. In the case of N branes, we take W = 
W + g(T)(l + SJ^i^f)- To simplify the analysis, we 
consider examples where a, /?, 7 are chosen in such a 
way that the full potential continues to have a global 
minimum at ijj = ?po, where V(tpo) is assumed to vanish: 



SV ~ NS 



(y> 2 -^ 2 ) 2 

(2a - N^ 2 f 



(20) 



We have included the factors of N appropriate to the 
multibrane case, when all branes are coincident. 

The second kind of correction involves the Kahler 
potential, which for N branes was taken to be K = 
-31n(T+f- V\ |^i| 2 ) [l]|. This was always understood 
to be an approximate expression, good when ip 2 <C cr, 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 



X X 

FIG. 11: Rescaled brane-antibrane potential i1'2'2fi including 
Kahler and superpotential corrections. 

but subject to higher order corrections in ip^, for exam- 
ple \ipi\ ■ The denominator in V + SV would then be 
replaced by 

(2a - Ni> 2 ) 2 -> (2d - N(iP 2 + d^)) 2 (21) 

As shown in Appendix B, the Kahler corrections also 
give contributions to the numerator of the F-term poten- 
tial, whose Taylor series starts with a term of order ip . 
For small ip 1 we can absorb this leading term into the 7 
coefficient of (|19|) . The kinetic term IC a i4> a (f> b also gets 
0(dtp 2 ) corrections, which can be inferred from eq. (|B6|l . 
but these are negligible whenever the branes are slowly 
rolling. 

By the above reasoning, the two kinds of corrections 
enter into the rescaled potential (JSJ) as 

( c( - 2 - 1)2+ ( 1+ (^fD 

(22) 

By experimenting with this generalized potential, one 
finds that the new coefficient c allows for the removal 
of the lower bound on b in l|10(l . Moreover, negative val- 
ues of d can remove the upper bound. It is easy to find 
values of c and d for which the potential has the desired 
behavior as a function of N, as in figure regardless of 
the value of b. For example c = 0.1, d = —0.5 yields the 
potentials shown in figure ITT1 at b = 0.01 and 6 = 5. It is 
important that the new parameters c and d are indepen- 
dent of the number of branes. Thus a continues to be the 
only parameter which changes due to a brane tunneling 
event. 

As in the previous section, we scanned over parameters 
for which A varies, hence the number of branes, while 
keeping b, c, d fixed, for the two examples with b — 0.01 
and b = 5. Figure ^] shows the number of e- foldings 
versus number of branes. The general growth of N e ~ 
\J~W C is evident, as it was for the simple model of the 
previous section. 

Figure ITS1 shows the correlation of dn/d\nk with n s . 
It shows that the small-6 model tends to have a red spec- 
trum, with n s < 1, and negligible running. The large-6 
model also has n s < 1 for most realizations, but for rare 
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log 1() N c 

FIG. 12: Log of number of e-foldings versus log of the critical 
number of branes, for the examples illustrated in figure 9. 




FIG. 13: dn/dlnk versus n s for the models with Kahler and 
superpotential corrections. 



cases with N < 20, it is possible to get large deviations 
from a flat spectrum with n s > 1. These are accompa- 
nied by sizable running, but as in figure ITU1 such points 
occur outside of the experimentally allowed region. The 
overall conclusion is that the generalized models have 
similar qualitative behavior to the simplest class which 
was studied in sections III-IV. 



VII. CAVEATS AND CONCLUSIONS 

In this work we have presented a very general means 
whereby a long period of inflation, with a suitably flat 
spectrum of density perturbations, can be automatically 
obtained after tunneling of branes, from a metastable po- 
sition in a Calabi-Yau manifold, into a warped throat 
containing antibranes. It depends only on the fact 
that the brane-antibrane potential generically has such a 
metastable minimum when the number of branes is suf- 
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ficiently large, and the possibility of starting at a place 
in the string theory landscape where the number of mo- 
bile branes exceeds the threshold for having a minimum. 
Once a certain critical number of branes N c is reached, 
the minimum disappears and the branes start rolling 
slowly to the true minimum. They automatically start 
from the same position, which is guaranteed to be close 
to the flattened region of the potential. 

In this scenario, an arbitrarily long time will be spent 
in a metastable de Sitter phase resembling old inflation, 
before the slow rolling period (similarly to ref. |l9j). Each 
tunneling event results in a small bubble of open uni- 
verse, within which another period of quasi-de Sitter ex- 
pansion will occur. The process repeats itself until the 
critical number of branes is reached. For us, these long 
periods of old inflation are irrelevant — we have focused 
entirely on the final phase of conventional inflation asso- 
ciated with the smooth motion of the branes until they 
annihilate with antibranes in the throat, hopefully lead- 
ing to the reheating of the universe 12(1 , and perhaps the 
production of cosmic superstrings |2l| . 

We numerically scanned over large ranges of the loga- 
rithm of relevant combinations of parameters which gov- 
ern the shape of the inflaton potential, and found that 
large fractions (~ 10 — 90%) of the trials resulted in in- 
flation with no need for tuning, whenever the flattening 
mechanism could occur. This seems to be a quite promis- 
ing approach for ameliorating the need for fine-tu ning 
within the framework of the string-theory landscape [23j . 

For convenience of numerical integration, we chose a 
form of the potential which has good behavior as the 
branes reach the bottom of the throat. This form, which 
differs from the commonly used Taylor-expanded form, 
is in fact the result of a direct perturbative compu- 
tation, and does not require any special assumptions. 
However, the perturbative derivation does assume that 
Ne/{ip — ipo) A <C 1 in order to be valid. One might there- 
fore wonder how sensitive our results are to the behav- 
ior of the potential when this approximation is break- 
ing down. We have investigated what happens if the 
Taylor expansion 1 — Ne/(ip — ipo) 4 is used in place of 
(1 + Ne/ (ip — ?/>o) 4 ) 1 , and we also checked different ways 
of regulating the singularity, like 1 — Ne/ (Ne+(ip — ipo) 4 ). 
For all these alternatives, we found that although the spe- 
cific parameter values where a flat potential occurs may 
be somewhat sensitive to the behavior of the potential 
near ipQ, the mechanism of multibrane inflation never- 
theless works in the same way as we have presented. We 
therefore believe it to be a robust result, which is qualita- 
tively quite different from field theory models of inflation. 

Although the dynamical flattening of the potential re- 
quires a rather special range of parameters in the sim- 
plest KLMT model, we showed that this restriction can 
be lifted once corrections to the superpotential and the 
Kahler potential are taken into account. We exhibited 
an example of such corrections that works, without sys- 
tematically exploring the parameter space of such cor- 
rections. Subjectively, it seems that the parameter space 



(referred to as {b, c} in section V) where multibrane in- 
flation works will be large, given that it was easy to 
find a successful example, but a more complete investiga- 
tion could be made. Moreover, such corrections should 
in principle be computed from a fully string theoretic 
starting point, which might restrict the allowed ranges of 
{b,c}. 

Another issue which deserves further study is the pre- 
cise relationship between the radial coordinate r of the 
brane within the warped throat, and the field ip which 
is used to describe the inflaton position in the super- 
gravity description. The exact correspondence between 
these fields was not needed in KLMT because the main 
dependence on ip was taken to be through the Kahler 
modulus factor (and superpotential corrections), rather 
than a combination of these with the brane-antibrane 
potential. We have simply identified the two fields in 
our analysis, but the actual relation could be more com- 
plicated, possibly introducing additional parameters into 
the Lagrangian. 
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APPENDIX A: BRANE-ANTIBRANE 
POTENTIAL IN WARPED THROAT 

We consider the background of AdS$ x A5 space with 
the AdS§ metric: 



R 2 



ds 5 = ^2 + d * 2 ) + -^2 dr 



(Al) 



where R 4 — Airag s Na' 2 is the radius of the AdS$ space 
created as the near-horizon geometry of a stack of N ^S> 1 
D3 branes. There is also an RR background field created 
by the stack of branes: 



A x o x i x 2 x 3 — lim 



1 



(A2) 



If one now writes the action for a probe brane is this 
background, one finds: 

S m - ~T 3 J d 4 ^- det (d^x M d e x N g NM ) 

+ T 3 Jd 4 ^ xW (A3) 

We will use the convention that £° = t, but the rest 
of the coordinates will describe the location of the brane, 
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including fluctuations (local bendings), x % = x l £M — 
x l (t, £ l J . We make here a further assumption, namely 

I 



that only its radial location varies. 



d^x M d{bX N gNM — d^aX°d^bX goo + d^x l d^,x l ga + d^rd^rg r 



The metric will look like: 
/ 1 



dr dr 
yrr q^q Qgi 



900 + 9rr ( Jjni) 

dr dr , ( dr \ dr dr 

yrr ~qFT ^TT hfxx ~~r y T r 1 g|T I yrr ~Q£T ~q£I 

9rr-§^-§gT Qyy + 9rr (jfiTZ^j 



dr dr \ 



dr dr 
yrrQ^lQ^ 

dr dr 
\ yrrQ^Q^) 



dr dr - > - 

'rr Qgo Q£2 yrr gga Q£3 

dr dr 
yrr-Q^TQ^ 

dr dr 

yrr&ri&rs 



g rr 



dr dr 

WW 



dr dr 

y rr WW 



9z 



g rr 



(*)' 



5oo 

fc o 

g yy 

g zz 



We can apply the usual formula: 



det (A + eB) = det A det (l + eA^B) ~ det A (l + e Tr (A~ 1 B)) 



2 

Writing the metric at the location of the brane as: g^ v = j^f^u we obtain: 



det (d^x M d^x N Qnm) = — 



R 



dr dr 
1 d^ 0£" 



(A4) 



dr dr 
"d^W 



(A5) 
(A6) 

(A7) 



r 



The factor in the bracket comes from the product 



(A8) 



Using these results in the DBI + CS action we obtain 
Sd3 - 



_ fA»f 



± T 3 / cT£ det 



9a; 



(A9) 



Expanding the square root term we obtain for a D3 brane 
(sign +): 



(A10) 



while for a D3 (sign — ): 



'D3 



T; 



2 J d^dZ" 



(All) 



The extra piece in the case of the D3 comes from the 
non-BPS nature of the antibrane, since the graviton ex- 
change and the RR exchange forces now add instead of 
canceling (for 3-branes the dilaton is constant and con- 
tributes no force). The potential energy for the D3 in 



the gravitational + RR field will make it favorable for 
the brane to move towards r = 0, that is, the bottom of 
the AdSs throat. For the setup in fjg.l |14|l we have to add 
the piece corresponding to the graviton + RR exchange 
between the D3 and the D3. However in this setup we 




r = 



Stack 




-D3 




FIG. 14: KLMT setup 



have to include the effect of the D3 in addition to that of 
the stack. The derivation of the background perturbed 
by the presence of the D3 placed away from the stack can 
be found in Appendix B of KLMT. It starts by writing 
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the metric in the form: 



by: 



ds 2 = ft-i (-dt 2 + dx 2 ) + hh [ dr 2 + ^g ab dy a dy l 



*- * 

~4 



(A12) 



placing the extra D3 brane at n and solving for the per- 
turbed harmonic function 



h(r) = ^+8h(r) 



(A13) 



which satisfies the Laplace equation in the space trans- 
verse to the brane, which is spanned by r and y a : 



which amounts to solving: 



I?' 



d r 



R 5 



d r 5h (r) 



= CS {6} (r-ri) 



(A14) 



(A15) 



which give 5h(r) ~ \r — r*i|~ 4 . A similar equation is 
satisfied by the unperturbed h. The constants C and C 
depend on the total "charge" , more precisely the total 
tension of the brane or stack of branes, so they are related 



c= c - 

N 



(A16) 



and since the solution for the unperturbed h (r) is already 
known, we can immediately find the solution for Sh (r) 



Sh (r) 



C 



r - r x 



N 



ril 4 



(A17) 



We can use the perturbed metric to calculate the action 
for the D3 using the new function h (r) + Sh (r) . 

Now modify the setup and place at ro a stack of M 
D3's with M <C N and an equal number of D3's placed 
at r\ , r2 . . . rj; , as shown in figure 1151 We can simply 
generalize the analysis above to come to the result: 



h(r) 



i? 4 R 



1 



M 

— - V 

N ^— ' \r ■ — r 



(A18) 



j 



There is no interaction between the D3's or between the 
Z33's in the second stack since they are mutually BPS. We 
can now go back to eq. (|A9|) and replace the expression 
i? 4 /r 4 by h(r). The relevant term is the second term 
from eq. IjAllj) . which in the limit 3> r becomes: 



2T 3 / d 4 ev-/ 



h(r) 



2T 3 f d^J^f^r r i S ft 



7? 4 



(A19) 



where ro is the location of the stack of antibranes. The 
above result is for one single D3 in the stack. The final 
form of the effective action should be multiplied by M, 
the number of antibranes in the stack. 




FIG. 15: Modified setup 



APPENDIX B: KAHLER POTENTIAL 
CORRECTIONS 

The Kahler potential is assumed to have the following 
form: 

K = -Mog(r) = -Mog(T + T-f(4>i^i)) (Bl) 

and we will also assume that / decomposes into a sum of 
functions, each one depending on one field <pf. 

K = Slog ^T + T-J^fii^i)^ (B2) 

The simplest form of the functions /j is 

h{<t>i,h)=M i (B3) 

but we can consider more complicated expressions, for 
example: 

fi(<(>i,0 i )=c(>i0 i -{<pi0 i ) 2 (B4) 
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Allowing for the more complicated expression for In the general case the Kahler metric derived from the 
fi [4>i,<f>A changes the structure of the singularities of Kahler potential of eq. I|B2[1 has the form: 
the Kahler potential and of the F-term. 



/ 1 — ^ 

/ 1 90i 

ah ah dfx , 



3fn- 



aK-i 

\ 90* 



au-i ah 

90*„! 90i 

ajn ah 

90* 90i 



ah 

d<f>2 
dfi dh 



g/a-i 9/ 2 

d<t> 2 
ajn 0f 2 

90* 902 



9/„ 
90„ 
9/i 9/„ 



9/^-i 9/„ 
d<f>* l _ 1 0(pn 
dfn df n , 9 2 /„ 
90* 90„ ^ 90„90*, 



(B5) 



The inverse can also be calculated: 



/C ah = - 



9/i dh I 9 2 /i , 
90i 90* / 90i90* "T" " "f" 
9/i / 9 2 /i 
O0i ' d<f>id(f>l 
dh I d 2 f 2 



dfn Of, 



9 2 /„ __ 
90 n 90* 90 



0/2 / 

9<A, / 



9/i / 9 2 /i 
/ d<fiid4> 

1 / 90i90* 





9/: 



1/ 



9 2 / 2 



902 90| 



1/ 

' 90290| 



gA. / 9 2 /„ 
90* ' 90„90* 






9/„ 



9 2 /„ 



V- 



The F-term is calculated using the usual formula: 

Vp = (x^AJ^DbW^ - 3 |t¥| 2 ) (B7) 

where 

/C = -31og(r) (B8) 
The covariant derivatives are expressed as: 



D T W = 


8tW 4 


- J¥<9 T /C = 


= 9 T W^ - 


(B9) 


D^W = 






Oipi r 


(BIO) 


D T W = 


%F-( 


-TF%/C = 


= OjW - ~w 


(Bll) 


D^W = 






94>* i" 


(B12) 
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Using the above expression for the inverse Kahler metric one obtains the following terms: 



K 1 1 D T WD T W = - d~W --W) [ d T W - -W 



3 V 



dfn dfn 



RT^DtWD^W 

k^ t d^Wd t w 

K^DZWD 6 W 



d-W --w\W 




r -(d T w- 3 -w)w i:u>L 



dfi q a , dfn 

o 2 h + r ; 



1 001901 




(B13) 
(B14) 
(B15) 
(B16) 



Adding all the above terms and —3 \W\ we obtain the F-term: 



V F = 



dfi dfi 

^WdrW [ 1 + 



+ ••• + 



dfn dfn 

d<j> n <90* 
r d 2 f~ 

d<t>nd4>' n 



- r (WchjW + Wd T W) 



(B17) 



We assumed that the superpotential itself is independent of the fields <f>i . . . t 
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